A DISSIPATIVE MODEL FOR HYDROGEN STORAGE: 
EXISTENCE AND REGULARITY RESULTS 



ELISABETTA CHIODAROLI 



Abstract. We prove global existence of a solution to an initial and boundary value problem for a 
highly nonlinear PDE system. The problem arises from a thermomechanical dissipative model describing 
hydrogen storage by use of metal hydrides. In order to treat the model from an analytical point of 
view, we formulate it as a phase transition phenomenon thanks to the introduction of a suitable phase 
variable. Continuum mechanics laws lead to an evolutionary problem involving three state variables: 
the temperature, the phase parameter and the pressure. The problem thus consists of three coupled 
partial differential equations combined with initial and boundary conditions. Existence and regularity 
of the solutions are here investigated by means of a time discretization-a priori estimates-passage to the 
limit procedure joined with compactness and monotonicity arguments. 



1. Introduction 

The paper deals with a thermo-mechanical model describing hydrogen storage in terms of metal hy- 
drides. Hydrogen storage basically implies a reduction in the enormous volume of hydrogen gas; 1 kg 
of hydrogen at ambient temperature and atmospheric pressure has a volume of 11 m^. There are basi- 
cally six methods in order to store hydrogen reversibly with a high volumetric and gravimetric density: 
hydrogen gas, liquid hydrogen, physisorption, complex, chemical hydrides, metal hydrides. This latter 
technique exploits the possibility of many metals to absorb hydrogen: such metals and alloys are able to 
react spontaneously with hydrogen and they can store a large amount of it. These materials, either a 
defined compound or a solid solution, are designed as metal hydrides: their use in the hydrogen storage 
is of interest in terms of safety, global yield and long-time storage. Energetic and industrial applications 
warrant this interest. Indeed, hydrogen is foreseen to be a clean and efhcient energy carrier for the future 
(cf. [TC]). Nowadays, energetic needs are mainly covered by fossil energies leading to pollutant emissions 
mostly responsible for global warming. Hydrogen stands among the best solutions to the shortage of 
fossil energies and to the greenhouse effect, in particular for energy transportation. 

Our research moves in the direction of providing a predictive theory for the storage of hydrogen by use 
of metal hydrides (cf. [H]). To this aim our analysis refers to the thermo-mechanical model introduced 
by E. Bonetti, M. Fremond and C. Lexcellent in [7] and it complements their results. Following the usual 
approach of Thermodynamics, Bonetti, Fremond and Lexcellent have derived the governing equilibrium 
equations. The analytical formulation of the problem they obtained was new and, to the best of our 
knowledge, no other related results can be found in the literature. Our work is attempted to go one step 
further both in modeling and analytical aspects of hydrogen storage by use of metal hydrides (cf. [5]). 
The problem we will deal with results from a phase transition model and is formulated via the classical 
principles of Continuum Mechanics. The related strongly nonlinear PDE system has been investigated, 
from the point of view of existence and regularity of the solutions. 

In order to get acquainted with the phenomenon, we recall some of its basic features. Some metals are 
able to absorb hydrogen atoms and combine with them to form solid solutions. We assume the existence 
of two solid solutions: the a-phase and the /3-phase. The presence of one phase with respect to the other 
depends on the pressure of the hydrogen. To provide a good mechanical model to be analytically treated 
from the point of view of existence and regularity of solutions, it seems useful to exploit the theory of 
phase transitions. We choose the volume fraction of one of the phases as a state quantity and we denote 
it by X- Hence, x satisfies the relation 

XG[0,1], (1.1) 
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and, assuming that no voids appear in the mixture, the volume fraction of the other phase is simply given 
by 1 — x- More precisely, if x = 1 we have the a-phase, if x = we have the /3-phase, and if x € (0, 1) both 
phases are present in suitable proportions. The state variables of the model are the absolute temperature 
9, the hydrogen pressure p, and the phase parameter x along with its gradient Vx accounting for local 
interactions between the different phases. Constitutive relations for the state quantities will be chosen 
in such a way that the principles of Thermodynamics are satisfied. Finally, to describe the thermo- 
mechanical evolution of the system, a pseudo-potential of dissipation will be considered. The constitutive 
equations will be substantially recovered as in [7]. Nevertheless, some improvements will be introduced 
in the model. Indeed, in [7] the PDE system is written by neglecting dissipative effects and microscopic 
movements in the power of interior forces. Our approach will be different, as we aim to derive the complete 
dissipative model accounting for microscopic velocities and diffusive phase transformations. Moreover, 
[7] is concerned with the study of a weaker formulation of the system: in this framework the authors have 
been able to prove a global existence result holding for n = 3. We complement their results. Indeed, in 
our framework we can prove a more general existence result in the three-dimensional setting. Moreover, 
by refining the assumptions on the data, we are able to show further regularities for the solutions to our 
problem. Some related analytical results can be found for models describing irreversible phase transition 
phenomena. In particular, we refer to [T^ and [T^: these papers deal with nonlinear systems of PDEs 
governing the evolution of two unknown fields (0 and x) 3'iid prove some existence results. We refer also 
to [5] and [6], both concerning the analysis of a dissipative Fremond model for shape memory alloys. 
The problem therein is quite similar to our both for modeling aspects and analytical investigation: many 
techniques and analytical tools from [5] and [6] have inspired our proofs. 

Now, let us introduce the complete thermo-niechanical model describing hydrogen storage by use of 
metal hydrides and including dissipative effects and miscroscopic velocities in the constitutive equations 
as well as microscopic forces in the principle of virtual power. Then, we state an initial and boundary 
value problem for the obtained model. 

We write the model in terms of the state variables 9, p and x- In order to achieve a precise description 
of the phenomenon, we take into account both the equilibrium and the evolution of the system, which are 
characterized by two functionals: the free energy ^E*, defined on state variables, and the pseudo-potential of 
dissipation $, defined on dissipative variables. Constitutive relations for the involved thermo-mechanical 
quantities are chosen in accordance with the principles of Thermodynamics. 

We preface our discussion by defining the hydrogen density pH and the total density p (in what follows 
we take p = 1). Thus, letting ch be the capacity of the hydrogen, we have that ch = Ph{p — Ph)~^, i-e., 

CH 1 

PH = P— = -, 

I + Ch t 

where t stands for the hydrogen's specific volume. 

Now, we introduce the free energy function '^'{9, r, x, Vx)- By thermodynamical and duality arguments 
it follows that the free energy is concave with respect to 0, while we assume that it is convex with respect 
to T, X, and Vx- Thus, as ^ is convex with respect to r, we can introduce the dual function of 4* as 
follows: 

**(^?, C, X, Vx) = sup{Cs - ^{0, s, X, Vx)}. 
Next, we define the pressure p connected with the hydrogen's specific volume t by the following relation 

a* 

-p:= — (0,T,x,Vx). (1.2) 



This corresponds to setting 



(0,-p,X,Vx) --^(f?,-p,X,Vx). 



d{—p) ' ' ' Qp 

Thus, if we assume sufficient regularity for the functionals we get 

4'*(6',-p,x,Vx) = -pr-*(6l,r,x,Vx). 

Finally, we deal with the enthalpy functional G{9,p,x,^x) defined in terms of the Legendre-Fenchel 
transformation of with respect to the specific volume r. More precisely, the enthalpy G{p) is defined 
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by 

G(p) := -v&*(-p) 

so that it results 

G(0,p,X,Vx) - ^'(0,T,x,Vx) +PT. (1.3) 
In particular, we recover that G is concave with respect to p and 9, while it is convex with respect to x 
and Vx- 

Then, we make precise the constitutive relations holding for the entropy s, the specific volume r, and 
the internal energy e (see (|1.3p l. We have 

a* dG ,^ ^, 

dG , , 

e^'$ + es = G-pT + es. (1.6) 

The fundamental balance laws of Continuum Mechanics, written in a smooth bounded domain Q C K" 
during a finite time interval [0, T] are: the momentum balance, the energy balance and the mass balance. 
In the following we will denote the time derivative of any function / by the symbol ft. At first, by the 
principle of virtual power written for microscopic movements (i.e., neglecting any virtual macroscopic 
velocity), we recover an equilibrium equation for the interior forces, which formally correspond to the 
balance of the momentum seen as a microscopic equilibrium equation. Namely, in absence of external 
actions, we readily get 

B-divH = infl, (1.7) 

coupled with the boundary condition 

H n^O on r, (1.8) 

where F = 50. From (|1.7p - (|1.8p we can deduce the mechanical meaning of the vector H, which indicates 
a work flux vector, while i? is a scalar quantity collecting microscopic forces. In case when macroscopic 
deformations are described by — rt, we address the following energy balance equation 

et + divq^r + Bxt + H -Vxt-pTf in f7, (1.9) 

where by e we denote the internal energy as defined in (jl.6p : q represents the heat flux vector, for which 
we will state later a suitable boundary condition. The right hand side of (|1.9p accounts for heat sources 
induced by mechanical and external actions. More precisely, r is an exterior heat source, while heat 
sources induced by microscopic forces are collected by the terms involving the quantities B and H. The 
presence of microscopic mechanically induced heat sources in ()1.9p is justified by a generalization of the 
principle of virtual power in which interior microscopic forces and motions are also considered, as they are 
responsible for the phase transition (cf. [12]). Finally (|1.9p is complemented with a Neumann boundary 
condition 

-q-n^O onF, (1.10) 

where n stands for the normal unit vector on the boundary F = dil. This corresponds to prescribe 
a null heat flux through the boundary F. Here, and in the remainder of the work, we assume small 
perturbations. Hence, letting the mass of the hydrogen that is not in the solid solutions keep constant, 
the hydrogen mass balance reads as follows 

{pH)t +divv = infl, (1.11) 

where v is the hydrogen mass flux. Then we combine (|l.lip with the following boundary condition 

-v-n + jp = onF, 7 > 0, (1.12) 

by which we require that the hydrogen flux through the boundary is proportional to the difference between 
the exterior and the interior pressure (here the exterior pressure is chosen equal to 0). 

In order to describe the thermo-mechanical evolution of the system and to include dissipation in the 
model by following the approach by Moreau (cf. [20]), we introduce a pseudo-potential of dissipation $ 
depending on xt, Vxt and W9. We recall the properties of <&: 

$ > 0, $(0) = 0, $ is convex. (1.13) 
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By (|1.13p . it turns out that the subdifferential 9$ is a maximal monotone operator with € 9$(0). In 
particular, it follows that 

a<f>(xt, vxt, V0) • (xt, vxt, , ye) > o. (i.u) 

Now, we are in the position to exhibit the constitutive relations for B and H, in terms of G and $. 
Unlike Bonetti, Fremond and Lexcellent (see [7]) we include dissipative effects both in B and in H and 
we prescribe them to be given by the sum of a dissipative and a non-dissipative contribution. Let us 
introduce a useful notation: nd in the apex is used for pointing out nondissipative contributions, while d 
stands for dissipative ones. Namely, we specify B as 

S = B-' + B^' = |^ + |^, (1.15) 

and H as 

The heat flux vector q is assumed to fulfil the standard Fourier law 

q = -koV0, (1.17) 

where ko > 0. Let us anticipate that, by a suitable choice of $, the heat flux can be expressed by use of 
the pseudo-potential of dissipation. Finally, we set the following relation for the hydrogen mass flux 

V = -AVp, (1.18) 

for A > (take, e.g., A = 1). 

Let us come to the functionals that describe the equilibrium and the evolution of the system. As first 
we set the enthalpy functional G as 

G{0, p, x,r)^a logp + bx{\ogp - logp,) - c^O log 6 + ^ | Vx|' + /[oa] (x), (1-19) 

where Cp > 0, (5 > and /[o.i](x) := if x G [0,1] and /[o,i](x) •= +oo otherwise. By p^ in (|l.f 9p we 
denote the equilibrium or Plateau pressure which is strongly temperature dependent. In accordance with 
physical experience, we let a > and b > (take, e.g., a = 6 = 1). Hence, experiments show that for 9 
sufficiently large the Van't Hoff law holds (see [12]), i.e., 

logpe = -cii + C2, (1.20) 

where Ci and C2 are positive constants. However, as we have already pointed out, the enthalpy G has to 
be concave with respect to the temperature, on the whole temperature interval. Thus we complement 
relation (|1.20p by setting 

logpe = hie), (1.21) 

where /i is a sufficiently smooth function, e.g. 

h{9) = —ci9^^ + C2 for 9 sufficiently large, say 9 > 9^, 
h{9) = h(9) for 9^^<9<9^, 

h{9) = constant for 9 < 9„, (1.22) 

with h{9) suitably defined to yield h e C^(R). Moreover, h{9) and other values have to be chosen in such 
a way that 

^ = -bxh"{e) - I < 0. (1.23) 

Indeed, (|1.23p ensures that G is concave with respect to temperature, which is necessary to get the 
thermodynamic consistency of the model. As it will be clear in what follows, (|1.23p is in direct relationship 
with some assumptions on h concerning the analytical solvability of the resulting heat equation. 

In order to derive the model in a thermodynamic frame we have to fulfil the second law of Thermody- 
namics in the form of the Clausius-Duhem inequality. To this aim, it is convenient to introduce the heat 
flux vector q formally as a dissipative quantity defined by the pseudopotential of dissipation $. Thus, 
we set 

1>=|lx*r + ^|Vx*|V||V^r (1.24) 
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for fco > 0. Hence, we define a new dissipative quantity 
so that letting 

q = 6Q'^ (1.26) 
yields the Fourier relation (|1.17p . Now, we point out that by use of the chain rule in (|1.9p and the 
above constitutive relations (|1.4p - p.6p . p.lSp - p.lSp . (|1.25p - (|1.26p . we can equivalently rewrite the energy 
balance (|1.9p as 

e (st + div Q'' - ^) = Q"' ■ve + B'' ■xt + H'' ■ Vxt- (1.27) 

Thus, after observing that by (|1.15p . (|1.16p and p.25p the right hand side of p.27p corresponds to 

9$ (9$ (9$ \ 

we get, thanks to (|1.14p . 

-Q''-V0 + S^-Xt + -H''^-Vxt>O. (1.29) 
Finally, dividing (|1.27p by the absolute temperature 9 yields the Clausius-Duhem inequality ensuring 
thermodynamic consistency, namely 

St + divQ'^-^ >0. (1.30) 
u 

Remark 1.1. Indeed, we have to remark that all the addends in (jl.29p turn out to he non-negative. In 
particular, we have 

-Q'^-Ve>0 (1.31) 

B^-Xt>0, (1.32) 

H''-yxt>0. (1.33) 

Thanks to (jl.3ip - p.32p . if we take v ~ Q in p.24p . the model still fulfils the Clausius-Duhem inequality. 
Indeed, letting v = Q in (|1.24p corresponds to neglecting dissipative effects in H (see (|1.16p ); this choice 
would lead to the model derived by Bonetti, Fremond and Lexcellent in [7] which actually is thermody- 
namically consistent. 



Now, we can write the system of PDE's in terms of the unknowns, by substituting in p.7p . (|1.9p and 
(jl.lip . the constitutive equations of the model specified by enthalpy and pseudo-potential of dissipation. 
In particular, we write the vectors of microscopic forces (jl.isp and p. 161) on account of (|1.19p and (|1.24p 

B = b(logp - h{0)) + a/[o,i] (x) + MXt, (1.34) 

if = (5Vx + i^Vxt. (1.35) 

Let us note that 9/[o.i] in (|1.34p stands for the subdifferential of the indicator function of the convex 
[0, 1] C K and it is obtained as a generalized derivative with respect to x of the non smooth function /[o,i] 
in (|1.19p . Actually, in our analysis we will consider a more general maximal monotone graph (3 in place 
of 9/[o,i], still with non negative and bounded domain. Thus, the complete PDF system originating from 
prg]) . fT7TT|) and JLT]) is written in Q := D, x (0,r) as follows: 

ibh"i9)0x + cp)0t - koA9 = r + nxl + y |Vxt|' - h0h'{e)xu (1.36) 

liXt - v^Xt - <5Ax + C = -h{-h{6) + logp), (1.37) 

i e /3(x), (1.38) 
V 

^a + bXy 

Then (|1.36p - (|1.39p are combined with initial conditions 

0{O) - 00, (1.40) 

X(0)=XO, (1.41) 

p{0)=Po, (1.42) 



XAp = 0. (1.39) 
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and the natural boundary conditions ((?„ is the normal derivative operator on the boundary T) 



kodnO = 0, 

{dnX)t + ^9nX = 0> 
XdnP + 7P = 0, 



(1.43) 
(1.44) 
(1.45) 



on r X (0,r). Observe that the unusual boundary condition (|1.44p results from the position (|1.35p and, 
however, it is perfectly equivalent to the standard boundary condition dnX = whenever the compatibility 
condition dnXo = holds true on F, as in our approach (cf. (12.34^ 1. 

Let us point out to the reader that we deal with a slightly modified version of the system (|1.36|) - (|1.39|) . 
as in equation (|1.36p we consider i/ = and neglect exterior heat sources (r = 0). From a mechanical 
perspective, letting v = in (|1.36p corresponds to require that dissipative effects on the gradients of the 
phases are negligible in the power of interior forces with respect to the other mechanically induced heat 
sources, which is reasonable in the framework of the small perturbations assumption. Thus in the sequel, 
by abuse of notation, we will refer to (|1.36p but considering v = Q and r = 0. Moreover, as our existence 
theorem for weak solutions to problem (|1.36p - (|1.45p does not ensure the pointwise validity of (|1.37|) . we 
have to consider an extension of inclusion (|1.38|) in the framework of a pair of spaces in duality. 

The above formulation of the problem of hydrogen storage complements the model advanced in [7]. 
In particular, the quadratic term in ()1.36p is new, as well as Axt in (|1.37p . Let us note that the 
solvability of the resulting PDE system, written as a phase-field problem, turns out to be interesting 
from the analytical point of view. Indeed, the system (|1.36p - p.45p is highly nonlinear and solving it (in 
some suitable sense) requires non trivial analytical tools. This is mainly due to the coupling of higher-order 
nonlinear contributions involving the unknowns, a maximal monotone graph and a quadratic dissipative 
term for the phase parameter, a logarithmic term involving the pressure. More precisely, in the parabolic 
equation (|1.36p the specific heat is a nonlinear function: to ensure coerciveness, we need to prescribe a 
suitable assumption on the function h. Dealing with the equation governing the evolution of the pressure, 
we have to combine the regularity of the function x and the pressure p, mainly to control the nonlinear 
evolution term. Finally, it is worth observing that the pressure has a major role in the evolution of the 
phase through the presence of its logarithm as a source in the corresponding evolution inclusion (see 
()1.37p ). The logarithm is easily controlled for high values of the pressure, whenever we are able to control 
p, but it degenerates as p \ 0. Thus, our proof present some ad hoc estimates for (|1.39p to handle this 
nonlinearity in the phase equation. 

Here is the outline of the paper. In the next section, we introduce an equivalent abstract formulation of 
the n-problem (|1.36p - (|1.39p . on account of the initial and boundary conditions (|1.40p - ()1.45p . In particular, 
we will focus our investigation on the three-dimensional problem as it is more meaningful from the physical 
point of view. At first, under suitable assumptions on the function h in (jl.2ip . we can state a global 
existence result fTheorem l2.3p holding for = in (|1.36p . Finally, under refined assumptions, we establish 
some further regularities for the state variables and the positivity of the temperature (Theorems 12. 4[ 12.51 
andES]). 

The existence result is proved in Sections 3 4 - 5 by exploiting a semi-implicit time discretization 
scheme combined with an a priori estimate-passage to the limit procedure. In section 4, we prove the 
positivity of the temperature and further regularities for the state variables, by performing suitable a 
priori estimates on the solutions of the problem. In the Appendix, we present some results we will refer 
to in our investigation. 



Our analysis refers to an abstract version of the problem (|1.36p - (|1.45p . We render the physical constants 
to 1 (i.e., a = 6 = Cp = A = fco = '5 = M = 7 = l)- Next, we start by listing the main mathematical 
hypotheses of our work. At first, we want to specify (3: we introduce a convex, lower semicontinuous and 
proper function 



2. The abstract problem: main results 



/3 : R (-00,-1-00; 



(2.1) 



satisfying the following property 



iniD{P) = (0,A^) 



(2.2) 
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for some constant > 0. By D{/3) we denote the effective domain of /?. In addition, we assume that 
there exists x* € (0, A^) such that 

= ${x*)<m Vsei?(/3). (2.3) 

Now, we can introduce the subdifferential of /3, i.e., 

/3 = 9/3; (2.4) 

then, /3 turns out to be a maximal monotone graph in M x R (we refer to [3] and [8] for basic definitions 
and properties of maximal monotone operators). Let us note that (j2.3p impHes € /3(x*)- In the above 
positions D = D{I3) denotes the domain of the graph /3, i.e., D = {r G R : /3(r) ^ 0}; we recaU that 
D g D0). 

Remark 2.1. The graph /3 is introduced in order to generalize the graph i9/[o^i] in the physical model we 
derived; in particular let us stress that the choice jS = 9/[o.i] would force x to attain values only in [0, 1] 
(see (|l.ip ]■ Nonetheless, more general classes of constraints could be admissible for the phase variable 
than that provided by ^ £ c^-^[o.i] (x)- particular, from a mathematical point of view, in this analysis 
we are able and aim to deal with an arbitrary maximal monotone graph (5 with domain included in some 
bounded interval [0, A^]. 

According to Remark 12.11 in (|1.38p we have considered /3(x) in place of 9/[o,i](x). Thus, the system 
p.36p - p.39p can be rewritten in Q as follows (recall that v = Q and r = in p.36p ). 

{h"{9)ex + m ^Ae^xt- oh'{9)xt, (2.5) 

Xt-iyAxt-Ax + ^ = hi9)-\ogp, (2.6) 
?€/3(x), (2.7) 
_P_ 

For the sake of clarity, before proceeding, we introduce some useful notation. We set 

and identify H with its dual space H' , so that 

V ^ H (2.9) 

with dense and compact injections. Let (•, •) and |j-|j be the inner product and the corresponding norm in 
H, and denote by (•, •) the duality pairing between V' and V. Hence, we introduce the following abstract 
operators: 

A:V^V', {Av,u) ^ [ \/v -Vu, u,veV (2.10) 



-Ap = 0. (2., 



Jn 

B-.V^V, {Bv,u) = [ VvVu+ [ vu, u,veV. (2.11) 

Then, to simplify notation, we set 

W := {/ € H^{n) : d„ f = on T}. 

Clearly, we have that W C V with compact embedding. Whereas the results we present refer also to 
lower dimensional space domains we have restricted our investigation to the three-dimensional framework 
as it is more meaningful from the physical point of view. Thus, let 

rJcR^ and g:=Ox(0,r), (2.12) 

where T is a fixed final time. We assume 51 to be a smooth bounded domain. We associate the functional 
Jh and J on H and V to the function /3, as follows 

Jh{v) = / (3{v) ifveH and (3{v) G L^^l), (2.13) 



Jh(w) = +oo if V e i7 and /3(w) ^ Li(r2), (2.14) 
J{v)^Jh{v) ifweK (2.15) 
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As it is well known, and J are convex and lower semicontinuous on H and V respectively. Note that 
they are also proper since V contains all the constant functions. We denote by OhJh : — > 2^ and by 
dv,v'J : V 2^ the corresponding subdifferentials. We remind that 

C edHJHix) if and only if ^ e H, 

X&D{Jh), and Jh{x) < {Lx - + Jh{v) e H (2.16) 
C €dv,V'Jix) if and only if C e V , 

xeDiJ), and J{x)<{^,X-v)+J{v) V« G F (2.17) 

where D{-) denotes, as above, the effective domain for functionals and multivalued operators. Note that 
OhJh : H ^ 2" and dv,v'J : V 2^' are maximal monotone operators. Moreover, observe that for 
XA& H we have (see, e.g., Ex. 2.1.3, p. 21]) 

£, G dnJHix) if ^-nd o'^ly if ^ £ Pix) almost everywhere in O. (2.18) 
On the other hand, one can easily check that the inclusion 

^HJH{x)'^HC^^v.v'J{x) VxGF (2.19) 

holds, just as a consequence of (|2.15p (compare (|2.17p with (|2.16p '). Therefore, for x G and E H, the 
condition 

C G /3(x) almost everywhere in (2.20) 

implies 

C^dv.yJix)- (2.21) 
In addition, if we assume % G ^ G and (|2.2ip . then (|2.20p holds (cf. [3]). In particular, on account 
of (|2.19p . we achieve the validity of the following equality 

dHJHix)=Hndv,v'J{x) VxGK (2.22) 

Nonetheless, we should observe that (|2.22p is false for more general functionals, as one can verify by 
referring to the example presented in [?]. 

Now, we are in the position to rewrite the PDE system (|2.5p - (|2.8p . combined with initial and boundary 
conditions (|1.40p - (|1.45p . in the abstract setting of the triplet (V, H, V'). We have to remark that we are 
not able here to deal with a strong version of p.5p - ()2.8p . In particular we cannot deal with the natural 
extension of /3, i.e. the subdifferential OhJe- Hence, the variational inclusion governing the dynamics of 
the phase is written in the abstract setting of the V' — V duality pairing. Nonetheless, even if it cannot 
be stated a.e. in Q, it retains its physical consistence since it forces the phase to attain only meaningful 
values. Indeed, if ^{t) G dv,v'J{xit)) fo^ almost any t G (0,T), we have in particular that 

x{t) G D0) a.e. in fl. (2.23) 

For instance, if we take /3 = -^[o.i]; tti^ abstract relation ^(t) G dv,v' J{x{t)) for s-^- t (0;^) implies 
that x{t) G [0, l]vQ, and consequently x G [0, 1] a.e. in Q. Thus the system is rewritten, in V' and a.e. 
in (0,r), as follows: 

et+Ae = -hie)xt+xl (2.24) 

e = 0-xih{e)-0h'{e)), (2.25) 

Xt+A{i^Xt+x)+^^h{e)-logp, (2.26) 

^edvyJix), (2.27) 

Ut+Bp = 0, (2.28) 

u^-^. (2.29) 
1 + X 

Let us note that, by (|2.25p and (|2.29p . we have introduced two new auxiliary variables: e and u. The 
variable e can be expressed as a function ip of the variables 6 and x- 

e = i>{e, x)-^o- x{h{e) - eh'{0)). (2.30) 



^By [0, l]v we denote the convex; [0, l]v := {v £ V : v ^ [0, 1] a.e. in Q}. 
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On the other hand, the variable u has a precise physical meaning: it can be interpreted as the normalized 
hydrogen density. Thanks to (|1.40|) - (|1.42p we are allowed to set the following initial conditions for e and 
u 

e(0) = eo := - xoiKeo) - Boh' {60)), (2.31) 
u(0) = uo := ^° ■ (2.32) 

Remark 2.2. We point out that, whenever ^ S dv.V'Jix) ^^'^ ^ H a.e. in (0, T) we have that actually 
^ G daJnix) '^•6- (0,T), from which one can deduce that ^ G /3(x) ^-e. in Q. 

Now, we set the hypothesis on the data prescribed in the first part of our analysis. Concerning the 
Cauchy data in (|1.40D - (|1.42p . we assume that 

6*0 e V, (2.33) 
XoeWnDiJH) (2.34) 
Po€V, logpoeL\n). (2.35) 

Note that (|2.35p yields po > Q a.e. in fl. Moreover, (|2.34p implies that xo G D{(3) a.e. in fl and in 
particular: 

< Xo < A/3 a.e. in (2.36) 
where the value is introduced in (|2.2p . Finally, from (|2.34p - (|2.36p we can deduce that 

uo e V. (2.37) 

In fact, the following estimate holds 

ll^^olly < C (iboll'^ + Ibolly llVxolly) < C, (2.38) 

where we have exploited the continuous embedding V C L^(ri). 

Hence, we ask for a suitable regularity of the thermal expansion coefficient h(9), in agreement with 
the assumptions leading to the physical consistence of the model (see (|1.23p ). We require 

hew^^°°{M.)nc^{R), (2.39) 
l|/i|liv^,=o(K) + |/i'(C)CI <c., \h"iCK\<cl VCeM, (2.40) 

for some positive constants Cfi,c'f^. In addition, let > such that (recall (|2.23p and (|2.40p ) 

1 + vh"{C)C >cs>0, V77 € D0), VC e M. (2.41) 

This correspond to assume that the product cJ^A/j is small with respect to 1. The hypotheses we made 
on h allow us to infer that ip (see p.30p ) is a bi-lipschitz function of the variable 9. Indeed the following 
bounds hold (see ([1301) -(EZD)) 

< Cs < dill; <1 + Ce, (2.42) 
\d2M < Ch, (2.43) 

for a positive constant Ce depending on cj^, where ditp, d2'>p denote the partial derivatives of ip with 
respect to first and second variable, respectively. Let us recall that for the initial datum eo we have: 
Co = "0(^01X0) (see (|2.3ip ). This and the regularity of 60 and xo hi (|2.33p - (|2.34p . along with the above 
properties of ip, easily yield 

eo e V. (2.44) 

By use of a semi-implicit time discretization scheme combined with the a priori estimate-passage to the 
limit procedure, we can prove the following related global existence result. 
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Theorem 2.3 (Existence). Let (EST]), (^Ili)) and ^^-^JQ hold. Then there exists 

a quintuple of functions {6,e,x^p,u) fulfilling 



9 e L°°{0,T,H)r\L^{0,T, V), 


(2.45) 


e G W^^\0, T, V') n i°°(0, T, H) n L2(0, T, F), 


(2.46) 


xeH\o,T, V)nL°"{0,T, W)nL°°{Q), 


(2.47) 




(2.48) 


p e i?i(o, T, ij) n i°°(o, T, V) n l2(o, t, iJ2(r2)), 


(2.49) 


logpeL°°iO,T,L\n))nL\0,T, V), 


(2.50) 


u e H\0,T,H)r\L^{0,T, V), 


(2.51) 



and solving (|2T24|) - ([2T29l) in V , a.e. in {0,T) along with (|1.4ip and ()2.3ip - p.32p . 

The regularities obtained in Theorem 12.31 allow the equation (|2.28p to make sense a.e. in il x (0,r): 
the same does not hold for (j2.24p and (j2.26p since the regularities in space are not strong enough. But, 
still dealing with the complete problem (|2.24p - (|2.29|) . we can establish some further regularity results for 
the state variables and, in addition, the positivity of the temperature 6. As first, we can prove some 
properties for the inverse of the pressure p. To this aim, we need to make an additional assumption on 
the inverse of the initial datum po- Namely, we need 1/po = Pq^ G H. Then, the following result holds. 

Theorem 2.4 (Regularity of the pressure). Under the same assumptions as in Theorem ] 2. SI if 

-=Po'^H, (2.52) 
Po 

then the following property holds 

p~^ e L°°{0,T,H)r\L'^{0,T,V). (2.53) 

The next result is concerned with improved regularities for the time derivatives of the state variables 
6 and x- 

Theorem 2.5 (Further Regularities). Let (j05)l - ([05)) . (EST]), (EH]), ((^3^ and ([1321) -(Elll) hold. 
Let the quintuple {9,e,x,P,u) fulfill (E^ - (E3T|) . (^35)) and solve ^^-^^ along with pnT|) and 

Xo e D{dv,v'J), (2.54) 

then it holds that 

e e H\0,T,H)r\L°°{0,T,V), (2.55) 

XeW^'°°{0,T,V). (2.56) 

Note that the improved regularity (|2.55p for the temperature 9 allow us to infer from (|2.25p that also 
the time derivative of the internal variable e belongs to L^{0,T, H). Therefore, the solutions obtained 
in Theorem 12.51 fulfil equation (|2.24p a.e. in x (0,T). Finally, we can establish the positivity of the 
temperature 9, holding under suitable assumptions on the logarithm of the initial datum 9o. Indeed, the 
following theorem holds. 

Theorem 2.6 (Positivity of the temperature). Under the same assumptions as in Theorem \2.5[ if 

\og9oeL\n), (2.57) 

then the following properties 

log^e L^{0,T,L\n))nL^{0,T,V) (2.58) 
^eL2(0,T,i/) (2.59) 



hold. 



Clearly, (|2.58p ensures the temperature 9 to be positive a.e. in Q. Theorems 12.41 [^75l and [^31 will be 
proved by performing some suitable a priori estimates on the solutions of the problem whose existence is 
stated by Theorem [ 
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3. Time discretization 



In order to achieve the proof of Theorem 12.31 we proceed as follows. First of all we establish a global 
existence and uniqueness result for an approximating problem. Next, we perform some a priori estimates 
that enable us to pass to the limit. 

In this section, we approximate the system (|2.24p - (|2.29p by use of a semi-implicit time discretization 
scheme. Thus, letting N be an arbitrary positive integer, we denote by r :— T/N the time step of our 
backward finite differences scheme. In the forthcoming analysis, we will extensively use the following 
notation. Let {V'^ , ...,V'^) be a vector. Then, we denote by Vt and Vr two functions defined on 
the time intervals [0,r] and (— oo,T] which interpolate the values of the vector piecewise linearly and 
backward constantly, respectively. That is, 

«.(0):-V^°, Vrit) := a,it)V' + {1 ^ a,{t))V'-^ (3.1) 

Vr{t) : = V° at <0, Vr{t):^V' ifte{{i-l)T,iT] (3.2) 

where 

a^(t): = li^ii^iMifte((i_i)^,^^], (3.3) 

T 

for i ~ 1, N . Moreover, let us introduce the backward translation operator TJ- related to the time step 
T by setting 

Trf{x, t) /(.x, t - t) for a.e. (x, t) e n x (0, T), 

V/ : O X {-T, T) ^ M measurable. (3.4) 
Next we regularize the initial datum for the internal variable p by defining 

J Po if Po > T- ,„ 
Pot = < (3.5) 
[ r lipo <T. 

Clearly, the regularized datum por is still in V. Let us note that, thanks to p.35p . logpor G L^(fi) and 
moreover 

borlly + ||logpor|lil(j2) < C, (3.6) 

for some constant C depending only on Upollyj l|logPo|lLi(j2)7 l^^l and T. Besides, in view of p.35p . (|3.5p 
ensures that 

logpor^H Vt (3.7) 

as well. Let us note in advance that, in the approximating form, we set the variational inclusion p.27p 
in H by substituting the abstract operator dy.v'J by the corresponding maximal monotone graph OhJh 
in H , provided we can prove some regularity of the solutions. As a consequence we will be able to solve 
the discrete variational inclusion a.e. in Q, (cf. also Remark l2.2p . Then, the approximated problem can 
be formulated as follows. 
Problem P-r- Find vectors 

(eo,e\....,e^) e (3.8) 

(x",x\....,X^)eW^^+\ (3.9) 
(FO,pi,....,P^) e (3.10) 

satisfying 

e°-0O, X"=X0, P"=Par, (3.11) 

and such that, by setting 

£;o = e°-x°(/i(e°) -e°/i'(e")), (3.12) 
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the following equations hold ior i = 1, N 





- E'-^ 




T 








- r-' 




T 







AQ' ^ -h{Q'-^)^ ^ +[^ <^ in V, (3.14) 



T 



(3.15) 



Ax' - Ax 



+ 1,^ ± +Ax"+^"--h{e'-^)-logP'-^ ill V, (3.16) 



T 



T 



BP' = in V', (3.17) 



for 



pi 

W = r, (3.18) 



E'edHJHixl- (3.19) 



Remark 3.1. Since we have set the variational inclusion (j3.19p in H and due to the regularity assump- 
tions, we will see that equation p.l6p turns out to make sense also in H . 

Let us note that p.l9p implies that S' G (3{x^) a.e. in il; in particular, since D{P) is included in some 
bounded interval [0, A^] (see (|2.2p and (|2.36p ). we can infer 

0<X'<A^ Vi = 0,...,7V, (3.20) 

a.e. in n. The conditions in ((XTT|) . combined with ([5TT^ and ^J^, provide that (cf. ((OT|l - ((0^ ) 

^° = eo, (3.21) 

U° = uor:=-^. (3.22) 

1 + Xo 

We can prove existence and uniqueness of a solution for the approximating discrete problem at any 
step T > 0. Indeed, the following lemma holds. 

Lemma 3.2 (Discrete well-posedness). Under the assumption (|2.33p - (|2.35p and (|2.39p - (|2.4ip . for 

any r > the problem Pr admits a unique solution. 

Proof. Owing to (|3.11l) - (|3.13p . ()2.33p - (|2.35p and (|3.5p . we can restrict ourselves to prove that for any 
fixed r > and for any i > 1, the system p.l4p - (j3.19p admits a unique solution. The main idea is to 
proceed by induction on i. Indeed, we suppose to know 

(0'"\x'"\^'"^) e X {Wr\D{JH)) X V. (3.23) 

with 

logP*"^ei/. (3.24) 

(see (12171)). We look for 

{Q\X\P')^V yi{Wr\D{JH))y~V (3.25) 
solving the resulting equations p.l4p - p.l9p and such that 

logP'eH. (3.26) 

We take first p.l6p and rewrite it as 

i i—1 

^ + fl + -) Ax' + = ^ + -Ax'-' + h{&'-') - \ogP'-\ (3.27) 

T \ T / T T 

Since we are assuming (|2.39p - (|2.40p and p.23p - p.24p we can observe that the right hand side, say is 
known in H. Thus relation p.27p can be equivalently rewritten as 

(r-iM + C + 9ff Jh)x' 9 J^, (3.28) 

where 

Cx' = (1 + >^T-')Ax\ 
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Actually, we would like to exploit the well known results on the maximality of sums of monotone operators, 
holding under particular regularity conditions. Namely, in this framework, we can invoke Theorem 17.21 
and get the required existence of a function G W O D{Jh) fulfilling p.l6p . On a second step, we take 
into consideration (|3.17p . where is now the solution of p.l6p . Equation p.lTp can be reformulated as 
follows 

pi pi — 1 

. + BP' = — (3.29) 

\ + r i + x'-^ 

The existence and uniqueness of a solution P"^ £V follows directly from the Lax-Milgram theorem, taking 
into account that x' satisfies p.20p . Before proceeding, we have to check that p.26p holds. To this aim, 
we refer to the First a priori estimate, where we will prove that log?/' £ V (see (|4.23p V In particular, 
this and the fact that obeys p.20p yield p.26p . Namely, owing to p.lSp . there is a constant C such 
that 

||logF*||^<C(||logC/*||^ + 11x^11 J. (3.30) 

Finally, letting in p.l4p x' be the unique solution of p.l6p and exploiting once more standard results 
on maximal monotone operators we can find a unique function 0* g solving the equation. Indeed we 
can rewrite p.l4p more explicitly in terms of 0' as 

T-^e* - T~^x'{K'd") - O'h'iQ')) + AQ' 

= T-^Q'-\l + h'{(d'-^)x''^) - T-^h{&-^)x' + (?^-^2^\ . (3.31) 



As the right hand side, say Q, is known in H the above relation can be equivalently rewritten as 

{T-\id + ni) + A)e' 3g, (3.32) 

where Id stands for the identity operator in H and 

7^,(e0 = x'(e^;^'(e0-/^(e^)) 

maps H into H. For the sake of clarity, in our notation the subscript i in 72., is used for pointing out 
the dependence of TZi on x''- Then, in view of assumptions (|2.40p - (|2.4ip it is not difficult to check that 
Id + TZi is a Lipschitz continuous and strongly monotone operator from H into H, and consequently 
coercive. Since ^ is a maximal monotone operator with domain W, the hypotheses of, e.g., O Cor 1.3, 
p. 48] hold and we get the required existence and uniqueness of a solution O* G to equation ()3.14p . 
This concludes our proof of Lemma 13. 2| since for any i, and any fixed r > 0, the corresponding triple 
{e\x\P') solves the system (|XTi)) - . □ 

4. A PRIORI ESTIMATES 

In this section, we aim to establish some a priori estimates on the time-discrete solutions whose 
existence has been proved in Lemma [3.21 By virtue of Lemma [3.21 and owing to the position p.ip - p.4p . 
we may introduce the piecewise constant and linear in time functions Or-, St, Xt, ^r, Pt, Ur, Or, e^, Xr, 
Pt, Ut interpolating the corresponding values. Thanks to these notations, the scheme p.l4p - p.l8p is 
restated as follows in V' and a.e. in (0, T) 

dter + A9r = -hiTrOr)dtXr + {dtXrf , (4.1) 

er=er-Xr{h{Or)-Orh'{er)), (4.2) 

dtXr + l^A{dtXr) + AXr + Cr = KTrOr) - Tr{\ogPr), (4.3) 

^redv.V'JiXr), (4.4) 

dtUr + Bpr = 0, (4.5) 

Ur - j^. (4.6) 

As concerns the inclusion (j3.16p we could write it in H in terms of the above introduced piecewise linear 
and constant functions. Nonetheless, in order to perform a passage to the limit procedure as r \, 0, we 
have to set this inclusion in the abstract framework of the V — V duality. Thus, instead of OhJh we have 
written the corresponding abstract operator dy.v'J in (|4.4p . Obviously, by the regularity of the solution 
vectors, the existence and uniqueness result we have proved in the previous section can be extended to 
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the abstract framework of V . In addition, let us observe that by construction (cf. p.ip ). Xt, ^t, Pt, e-r 
and Ur satisfy the natural Cauchy conditions (cf. p.Sp . (|3.11l) and p.2ip - (|3.22p ) 

Xr(0)=xo, (4.7) 
0r(O) = 00, (4.8) 
Pr(0)=POr, (4.9) 

er(0)=eo, (4.10) 

Ur{Q)^Uor. (4.11) 

Hence, we are going to prove some estimates on the approximating functions solving ()4.ip - (|4.6p : such 
estimates hold at least for t sufficiently small, but the involved constants do not depend on r. Indeed, 
our aim is passing to the limit in the above system as r \, 0, by compactness or direct proof, to get 
(|2.24p - (|2.29p solved in a suitable sense. Let us recall the trivial equality 

2a{a-h)=a^ + {a-hf -b^, Va, G R, (4.12) 

which will be applied in the following estimates on the discrete solutions. 

First a priori estimate. We first test equation (|3.16p by (x* — X*)i where the value x* is intro- 
duced in (|2.3p and satisfies e P{x*)- Let us observe that, by monotonicity of the operator ShJe, (|3.19p 
yields 



E\x'~X*)>0. (4.13) 

; 

Hence, exploiting the relation ()4.12p . by using (|4.13p . we get 

1 II 7 l|2 1 II i i-l l|2 1 II ,_i ||2 ly \ ||„ ,||2 

x'-x4H + T^h'-x' 'Wh-tt-Wx' '-x*L+ 1 + ^ II^X'L 



2^11^ ^"UH ■ 2t"'^ ^ iiff 2r ^11^ V 2t 



< / h{&-'){x' - X*) - I \ogP'-\x'-X*)- (4.14) 
Jn Jn 

Then, we would like to formally test equation p.l7p by 1 — l/[/*. Nevertheless, in order to make the 
desired estimate rigorous we perform a detailed procedure. Thus, let < e < 1 and 7e(-) be defined by 

iM := <! f (4.15) 




Next, we introduce the primitive function defined by 

/u 

7e(s)ds + l, (4.16) 

so that it results 

r£(u) = u — \ogu if It > e 

r.(u) = 1 + logi + ( 1 - i ) u ifM<£. (4.17) 



£ \ e ^ 

Now, we rewrite (|3.17p in terms of the variables ^-nd thus obtaining 



+ B{W{1 + X')) = Q- (4.18) 



T 

Hence, we test (|4.18p by 7e([/') and, thanks to the convexity of F^, we obtain 

- / T,{u^)-- f r,{W-')+ f \ww\'{i + xn%m 

T Jn T Jn Jn 



[ U'\/U'\/x"l'e{U")+ I U\l + x')le[U') <Q. (4.19) 
Jn Jt 
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-fT,iW)+f |VlogC/f + / W 

"T Jn Jnn{u^>e} Jrn{u^>a} 

-fr,{w-')+f {i + xl- f vlog[/*vx^ (4.20) 

Jn JrnW'>a\ JnnW^>e\ 



< 

T 



The second integral on the right hand side turns out to be uniformly bounded thanks to p.20p , while we 
can exploit the Young inequality in order to treat the last term of (|4.20p . Eventually, using once more 
(|3?20l) . from ([420| we can obtain 

- / r,(t/o + i / |viogc/f + / w 

Jn ^ Jnn{u^>e} Jrn{U'>o} 

<Co + - f r,([/'-i) + i / |Vxf . (4.21) 
T Jn ^ Jn 

where Cq depends on A^. Let us note that, by inductive hypothesis, we are assuming logP*~^ G H, which 
implies logJ/'"^ € H (see (|XT5)) 1. whence U'~'^ > a.e. in f]. Analogously, since P'-^ e V and ((X^ 
holds, we have also that U^~^ G H. In particular, we can infer that (U^~^ — logC/*^^) G L^{^1). Finally, 
we observe that the last integral on the right hand side of (|4.2ip is bounded thanks to the regularity of 
X* (x* € W). Thus, to pass to the limit as e \ in (|4.2ip . we can apply the monotone convergence 
theorem and get 



- I ([/'-log[/') + i I |VlogC/f + / 

<Co + - f [U'-' - logC/'-i) + i / |Vxf . (4.22) 
''' Jn ^ Jn 

Due to the previous estimates and to the Poincare-Wirtinger inequality, from (j4.22p we infer 

([/* - log G iSl) and log V G F, (4.23) 

whence [/* > a.e. in VL. Now. we combine the estimates (|4.14p and (|4.22p and sum them. By use of 
p.lSp we easily deduce 



1 II 2 ||2 1 1 i i 1 1 1 II % 1 11^ 

'1 V 



2 ' 2Ty " ^ 11^ 2t ^ 2t II 11^ 

+ - / ([/■'- log [7') + ^ / IVlogC/f + / V 
T Jn ^ Jn Jr 

+ I log(l + x'-^)x'<Co + i / (C/'-i-log[/'-i) 
Jn T Jn 

+ / M0*"')(X'-X*)- / \ogU'-\x' - X.) + I log(l+x'-')x*- (4.24) 
Jn Jn Jn 

Let us observe that the last integral on the left hand side is non- negative (see p.20p ). Besides, owing to 
p.20p , we can handle the last integral on the right hand side of (|4.24p as follows 

/ log(l + x'-^)x* < I X'~'x* < A| \n\ . (4.25) 
Jn Jn 
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Hence, by summing up (|4.24p for i = 1, m, with m < N = T/t, and multiplying by r, we obtain 
1 1 ™ 



2 

i=l 



+ f (C/™-logC/™) + iVT / |Vlog[/f +f]r / W 
< Co + i llx" - X*||L + ^ PX% + log t/") 

m „ m „ 

+ Vr / Me'-^)(X^-X*)- Vr / \ogW-\x'~X*), (4.26) 

where Co = (Co + A| Now, owing to (|2.39p and p.20p . by use of the Holder and Young inequalities 

we can find two positive constants Ci and C2 such that the following estimate holds 



1 1 ™ 

+ / ([/"-log[/'") + - Vr / |VlogC/'|' + Vr / 
<Ci + i||x°-x*||i + ^||Vx°ir^ 



+ Jf" + C2 J^^ |logf/°| + 2^/^ E ^ IIIoS^IIlHo) ■ (4.27) 

Let us remind that U'^ = uot = Pot/ (1 + Xo) ^^nd log U'^ ~ log uot = logpor — log(l + Xo)- Then, in virtue 
of (|2.36p and p.6p we can write 



1 1 .™. 

2 IIX" - X*IIh + 2 EHIVxIlJ, + ^ llVx"ll?^ 

+ / ([/-- log C/™) + -^r / |Vlog[/f + / C/^ 

Jn ^ rrr Jo ._i Jr 



m — l 



<Cs + l \\x" X4l + ^ IIVxlL + 2A, 5: r ||logC/^||,,(^) , (4.28) 

1=1 

for some positive constant C3 not depending on t. Moreover, let us note that the following inequality 
holds 

i([/* + |logJ7^|) log f/*). (4.29) 

As a consequence of (|4.29p and owing to p. lip we can finally obtain 

1 1 

2 llx" - x*llL + 2 E ^ P^fH + i 

1=1 

+ ^ \\Un\LHn) + I i|log^'"llLMn) + ^ E ^ /j V log f/f + E r ^ C/'^ 

-| m—l 

< C3 + - llxo - X*IIh + ^ llVxoIlL + 2A^ E ^ IHU'Wmn) ■ (4-30) 

1=1 

Thus, the Poincare inequality and the discrete Gronwall lemma [Ml Prop. 2.2.1] applied to (|4.30p imply 
(see ([2341) and (l236| ) 

+ II^TlrllLi(o,T,Li(r)) + l|l°g^TllL°°(o,T,Li(n))nL2(o,T,v) - C- (4-31) 
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In particular, note that m,- > almost everywhere in Q. 

Second a priori estimate. From the estimates on in (|4.3ip it is possible to deduce analogous 
regularities for — + Xt)- In particular, let us discuss the regularity of \ogpr- Since (|3.20p holds, 
we can infer that Xt > a.e. in Q and Xt € L°°{Q) (see (|2.23p and (|2.2p 'l. Hence, thanks to (|4.6p . we 
have that 

l|logPr|lioc(o,T,Li(0))nL2(o,T,y) 

< C4 ( ||log?2r|lLoo(o,T,Li(0))nL2(o,T,y) + \\Xt\\ L''°{0,T,L^{n))nL^{0,T,V} ) < C*- (4.32) 



Third a priori estimate. We test (jXTB)) by x' ~X'~^- Thanks to (|IT^ and (jXTg)) . we get 







2 

+ z/r 

H 


V(x' 




T 






V(x' - 




' 1 


2 


r 


H 2 



if 



9 II ^ \\h 



< [ {r-x'-'){h{e'-')-iogP^-'). 

Jn 



(4.33) 



Now, we can handle the integral on the right hand side of (|4.33p by means of the Young inequality and 
get 



/ {x'-x'-'me'-')-iogP^-') 

Jn 



< 



x' -x' ^ 


2 

H 




T 


H 


-5( 



\n\\\h\\ 



\\ogP' 



(4.34) 



Thus, we add for i = 1, ...,m, with m < N and we easily recover (cf. (P^ . (PTTj) and 



i-1 



1/ > T 

H 1=1 



if 



llNx-ll 



+ lj{x"l<C5 + l\\Vxo\\l+ I Pixo). 
Hence, owing to (|2.34p . we obtain the following bound 



\Xt\\h^O,T:V) + 



L=°(0,T;Li(n)) 



< c. 



(4.35) 
(4.36) 



Fourth a priori estimate. We proceed by formally testing p. 161) by t^x'. Due to the monotonicity of 
the operator OhJh, it turns out that r j^^ > 0. Hence, by exploiting once more (|4.12p and Young's 

inequality, by similarly proceeding as for (|4.33p . we owe to (|2.40p . (I4.3ip and write 



\ pX% - \ ||Vx'-ii + \ \Wfn \ \W-'tn + - \Ml 

< I \Ml + \ [\^\ + \HP'~'\\l) ■ (4-37) 

If we sum up in (|i37|) for i = 1, m, we get (see (P^ . (PTT|) and 

rn 

iivx'"!Ih + uril + E^ ll-^^'IlL ^ ii^^oIIh + PxoIIs + c,, (4.38) 

i=l 

for any m < N. Finally, the regularity assumptions (j2.34p on xo ensure that 

\\Xrh^^O,T,W)<C. (4.39) 
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Fifth a priori estimate. Test equation p.l4p by E'-. Exploiting once more (|4.12p and recalling 
that = x') (see !^JE\i ). we can write 

^ J- ^ II P« ?7^»-l|l^ llfi^*-l|P 



f |vef 9iV(e',xO+ [ ^Q'^x'd2tP{Q\xl 

Jn Jn 



= - / h{e'-')E' 



E' 



t X - X 



To handle the last integral on the right-hand side we use the extended Holder inequality and get 



E^ 



I X - X 



< E' 



X - X 



(4.40) 



(4.41) 



We recall the bounds on diip and 192-01 stated in (g^ll)- (11331) ■ Then, due to (H^OD, and the Young 

and Holder inequalities, we can infer 



1 II j||2 T 



E' - E 



H 



1 

27 



E^ 



I' +-llve^||' 



<C7||Vx'||^ + c,||i?^||^ 



x' - x' ^ 




\W\\h 


x' - x' ^ 


r 


H 




r 



(4.42) 



L-*(n) 



where Cy = c\/(2cs). Then, multiplying (|4.42D by r and summing up for i = 1, m, we get 





x'~x' ^ 




x' -x' ^ 








+ r 




^ ) 




r 


H 


T 





(4.43) 



In order to apply the discrete Gronwall lemma [141 Prop. 2.2.1], we have to treat the m-th term in the 
last summation of (|4.43p separately. Here, we sketch out such a procedure. Namely, we deal with the 
TO-th term by use of the Young inequality and get 

2 

X" - X"' " 

T 

H 



1 



L-*(0) 



(4.44) 



< - \\E"Xh + 2tc;; WdtXrWi.^^^TM) + 2 \\dtXr\\h(o,T,LHn)) ■ 

Finally, thanks to the Third estimate (see (|4.36p ') and to the continuous embedding V C L'*(r2), we 
conclude that there exists a positive constant Cg such that 

2 \ 



H 



-/'m ~/m— 1 

/V /V 



L*i(n) . 



<C78 + i||i?'""' 



H ■ 



(4.45) 



Owing to (jX^ . (10T|1 and (|05l) . 1^:1511 yields 

I IIE'-lli + %EHI™*Hi < a + - ||e„f„ 







2 

+ 


x' - x' ^ 


^ ) 


\^ 2 2 


T 


H 


T 





(4.46) 



Eventually, (|2.44p and (|4.3ip - (|4.36p allow us to make use of the discrete Gronwall lemma [HI Prop. 2.2.1], 
thus obtaining 

We recall that Ci- is defined as e^- = iP{Ot,Xt) (see (|4.2l) 'l. Thus, the relation between 6^- and Bt is 
bi-lipschitz (cf. (P^ - f^^ l and allows us to infer from that 

ll^^llL~(o,T,H)nL2(o,T,y) - (4.48) 
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Sixth a priori estimate. Here, we aim at deriving some fm-ther regularities for the variables 
and e-r from the previous estimates. At first, let us observe that (|4.48p . combined with (|2.42p - (|2.43p and 
with the regularity of Xt (see (|4.3ip 'l. yields 

l|er|lL2(o,T,y) (4.49) 

Indeed, since Cr = iP{0t,Xt), the following inequalities hold 

2 / II- ii2 



\\^r\\L^O,T,V) = J W^^tWh + J ||V6'r9lV-'(6'r,Xr) + VXrC'2V'(6'r,Xr)||^ 

< T ||e.||i,.(o,T,if) + Ce W^OrWl.^o.TM) + Xr\\ mo,T,H) < C, (4.50) 

where Ce depends on Ce (see (|2.42p V Moreover, from (|4.49p and (|2.44p we can also deduce that 



L2(0,T,V) 



< C. (4.51) 
Finally, let us observe that (|4.48p yields 

\\^^LHo,T,v')^^- (4-52) 
In addition, from (|4.36p we can easily deduce that 

Thus, by a comparison in (j4.ip . we can infer that 

\\9terhi^o,Ty') < C. (4.54) 



Seventh a priori estimate. We consider p.l7|) and we rewrite it in terms of P*, P* ^ and 
thus obtaining 



P' 



BP' = 0. 



(4.55) 



Now, we test (|4.55p by P'. Thanks to (|3.20p . there exists a positive constant = 1/(1 + A^) such that 

1 



1 + 



<1 \fi^O,...,N. 



(4.56) 



By standard algebraic calculations, we get 



1 



1 



2tJ^1 + X' 2t^1 + x*-i 2r 

,1-1 



(p.-i) 



1 



(P* - P* 



pi pi— I 



X - X 



1 



(l + X*)(l + X'-i) 2t 



1 + 



|VP'f 

,,i-l 



{PI 



i\2 



(i + xO(i + x'-')' 



Next, we observe that we can find a positive constant Cp such that 



Cp P' 



Hence, we multiply (|4.57p by t and get 



< 



VP' 



(pr 



(4.57) 



(4.58) 



1 
2 

< 



2 1 ^ (pj-l)^ --2 



1 + X' 



pipi — i_ 



n 

T~X 



1+X 

i-1 



T 

y 



l + X')(l + X'-i) 2 



p)i pi 



' 1 II -112 

-^+CpT PH' 

1 + X' II lly 



(4.59) 



To handle the first integral on the right-hand side of (|4.59p we exploit (|4.56p . the extended Holder 
inequality and Young's inequality to obtain 

2 



pi pi 



ii + r){i + r-') - 4 II 11^ 



A-l 



V 



|2 

\h ■ 



(4.60) 
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We proceed similarly for the last integral on the right hand side of (|4.59p so to get 



2 7o + + - 4 II ll^"^4Cp 



(4.61) 



Then, combining (|4.59p - (|4.6ip and smnming up for i = 1, ...m with m < N, we obtain (see p. Ill) ) 



E- 



pi pi 



H 



|2 

\v 



1 1 1 1 1 2 ^ \ — ^ 



F 



i=l 
^-1||2 



(4.62) 



Hence, we can apply the discrete Gronwall lemma [HI Prop. 2.2.1] and, due to p.6p . (|4.36p . deduce that 

l|PTllL°'(0,T,H)nL2(o,T,y) - (4.63) 



Eighth a priori estimate. We test (|4.55p by (P* — P' ^). By exploiting once more (I4.12p . Holder's 
and Young's inequalities and similarly proceeding as for (j4.49p . we write 



1 /■ (P* - P 



1 



-iiVPNI^ -i|IVP*"MI^ +i||p'll^ 
slr^llff 2 II llif+2ll ll^'(r) 



r.-l||2 



2 ir llL2(r) - ^ 



pi — l^pi _ 



j-1 



(i + xO(i + x'-')' 

The integral on the right hand side of (|4.64p can be estimated as follows (see (|4.56p ) 



(4.64) 



1 

T Jn 



< 



H 2C;^ 



X - X 



(4.65) 



By combining (|4.64p and (|4.65p , summing up for i ~ \, ...,m and exploiting (|3.1ip along with (j4.58p , we 
get 



1=1 
1 

2C. 



E- 



pi—1 


2 




T 


// 




x' - X 




2 


T 


V 



^P'^l<\\\^P,r\\l + \\\p,r\\l.^r) 



|2 



(4.66) 



By suitably applying the discrete Gronwall lemma [TH Prop. 2.2.1] it is not difficult to recover the 
following estimate (see p.6p ) 

(4.67) 



IIP 



lL°=(0,T,y) 



< c. 



Ninth a priori estimate. Now, we want to achieve the following estimate 

\Wr\\m(a,T,H)<C. (4.68) 

To this aim, it suffices to prove that H^r |1^2(o t //) — ^ ^'^'^ ll^t^^Tlli^i-g < C. Thus, in order to 
check the first estimate, we exploit (|2.37p . (|4 63p and have 

2 



^r\\L^{^0,T,H) 



Pr 



< 



WP^L'^iO.TM) - ^• 



(4.69) 



Now, let us deal with the second estimate. By virtue of p.lSp . (|4.56p . (|4.36p and (|4.67p . we can write 



||9tUi-||^2(o_T,//) - 2 {\\dtPT\\L2(^Q^TM} + II-P^IIl=^(0,T,V) \\'^tXT\\L^(o,T,V) 



< c. 



(4.70) 
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Thus, combining (|4.69p - (|4.70p we are able to deduce (|4.68p . In particular, since ||c?tUr|i/^2(-o t if) ^ C!, a. 
comparison in (|4.5p allows us to infer that ||'Bpr|1^2(-o t ff) — ^ ''^^'^ standard elliptic regularity results 
yield 

\\Pr\\LHo,T,HHn))<G- (4-71) 

In the end, the following estimate 

Wurh^^o.T.y) < C (4.72) 
is a consequence of (|4.36p and (|4.67p . since (see (|4.6p ) 



I'^tII L°°(0 TV) " 

(0,T) 



1 + Xr 



H 



PtVxt 



1+Xr {l+Xr? 



< C'lO \\Pt\\l^(o^t,V) (l + IIXr|ii«=(o,T, 



W) 



< C. 



(4.73) 



5. Convergence results 



In this section we aim to deduce some convergence results that allow us to pass to the limit as r tends 
to in (|4?T|) - (j46l) and ([47]) . (|4?T0|) - (|4lI|) so to conclude the proof of Theorem [231 Since we will obtain a 
great deal of convergences, we prefer to retrieve them step by step instead of presenting all the results in 
just one proposition. At first, we can combine the previous estimates ()4.3ip . ()4.32p . ()4.36p . ()4.39p . ()4.47p . 
gUl), ^Ml-^^, SIS, (mg), BTT|) . KT3 so to obtain, on account of (|3II1)-(E21), 



|Il°=(0,T,H') + \\Xr\\H^O,T,V) — 



-'■r\\L°^{O.T,H)nL^{0,T,V} 



-'T\\w^'^{0,T,V')nL^{0,T,V) 



< C 



r'^ W L-^ {O.T.H)nL'^{0,T,V) 
"-r|lL°°(0,T,y) + \\'^-r\\H^{0,TM) — ^ 

P-'-\\L'^{o,TX)nL'^{o,T,m{n}) + \\Pt\\h^o,t,h) — ^ 

logPr|lioc(o,T,Li(0))nL2(o,T,y) - 

for c not depending on r G (0, f), for a suitable f > 0. 

Now, it remains to pass to the limit in (|4?T|) - (|46| and ([4J| . (j4?T0|) - (|4?TT|) as t \ 
tation yields 

ilXr |Il'==(0,T,W) - llXoillv + llXr |Il«:(o,T,W) - 

and analogous estimates hold for ||er||ioo(-Q 7^ ||PT|lioo(Q y and ||ur||^oo(g 7^ y). 
well-known weak and weak star compactness results apply to ()5.ip - (|5.2p . (j5.4p - (|5.5p and ensure the 
following weak and weak star convergences to hold, possibly for a subsequence of r. 



(5.1) 

< c (5.2) 
(5.3) 
(5.4) 
(5.5) 
(5.6) 

0. An easy compu- 



(5.7) 



By virtue of ([SJj 



Xr^X inH\Q,T,V)nL°°{0,T,W), 

Cr^e ill L°°{0,T,H)r\L^{0,T,V), 

Ur-^u ill H\0,T,H)r\L^{0,T,V), 

Pr^P ill H\0,T,H)nL^{0,T,V). 



(5.8) 
(5.9) 
(5.10) 
(5.11) 



Let us stress that, even if we do not specify it, the convergence results have to be intended to hold up to 
the extraction of a suitable subsequence of r, still denoted by r for the sake of convenience. Then, owing 
to strong compactness theorems (see e.g. Cor. 4, p. 85]), by (|4.54p . (|5.8p - (|5.1ip . we get 



Xr-^X inC'{[0,nH^-'{n)), 
Cr^e inL^{0,T,H) 

Ur^u and Pr^P in C"([0,r],iJi-'^(17)), ife>0. 



(5.12) 
(5.13) 
(5.14) 
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Moreover, the following relation is fulfilled (cf. p.f p - (|3.2p ') 



\Xt - Xt\\l^(q^t,v)< max t 



2 

/ 

l<i<JV 



X 



2 

2 

(0 
H 



< T \\dtXr\\l2,Q,T.m ^ (5.f5) 



and analogous estimates holds for \\pt ^Pt\\'l^(oth) ^'^"-^ 11^^ ~ ''-'■■r\\'L°°{o t h)- While for the difference 
between and e,- we have 

\\St — erWj^l^Qj^y,^ < T II (9(61-11^1 < TC. (5.16) 

Finally, with the help of (|ET^ - (|?ll)l and ^J^-^J^, we are allowed to infer that 

Xr^X ini°°(0,r,y), xr^x inL°°(0,T,M^), (5.17) 
e^^e in L^(0,r, V^')> ^ e in L°° {0,T, H) n L^{0,T,V), (5.18) 

0^^6i in i°°(o,r,iJ)ni2(o,T,y), (5.19) 

Ur~^u mL°°{0,T,H), Ur ^ u in L°°{0,T,V), (5.20) 
pr-^P in L°°{0,T,H), pr ^ P in L°°{0,T,V) n L^{0,T, H^{n)), (5.21) 
logp, inL2(0,T,F), (5.22) 
where the last limit y will be identified in the sequel. 

Now, we want to improve the strong convergence for e,- in ()5.18p . To this aim, we first note that 

II^T ^ ^IIlp(o,t,v'') — W^'^ ~ ^IIl~(o,t,v) II^'t ~ ^IIli(o,t,v) ' (5.23) 

where, by virtue of (jS.lSp . ||e^ — ^I1l=°\o t V) — ^^^'^ II^'t ~ ^IIli(o t v) ^ ^ ^'^'^ all 1 < p < +oo. Thus, 
we have e,- — >■ e in L^(0, T, V'). In addition, we can perform the following estimate 

lls-r — e|li2(Q^y^^) < IIet — e||^2(Q^j^^y) ||ei- — e||^2(Q y y/) (5-24) 

and conclude that 

Br^e in i2(o,r,iJ). (5.25) 

Finally, since the relation between 9t and Bt is bi-lipschitz continuous (see ()2.42p ). from ()5.25p and (|5.17p 
we can deduce that 0t is a Cauchy sequence in L^{0, T, H) which is a complete space. Therefore, recalling 
the convergence Or ^ 9 in L°°{0, T, H) n L^{0, T, V), by uniqueness of the limit we conclude that 

Or ^9 inL^{0,T,H). (5.26) 

Now, by virtue of the above convergences ()5.25p . ()5.26p . (|5.17p and owing to the properties of ip (see 
(12321)- (ElSl)), it is easy to recover ([2:25| from gSl), i.e. e = ik{9,x)- 

Next, we deal with the logarithmic term in (|4.3p . Thanks to the strong convergence in (|5.2ip . there 
exists a subsequence still denoted by Pr such that 

Pt ^ P a.e. in Q, (5.27) 

and consequently such that 

logpi- logp a.e. in Q. (5.28) 

In principle, the a.e. limit of logpr could be — oo in a subset of positive measure (in which p = 0): but 
the property 

/ |logpr|^ < c Vr 
Jq 

and the Fatou lemma imply that 

|logp|^ < liminf / |logpr 
•^N-o Jq 

whence logp is well defined. Actually, we are now about to show that the weak limit y of the sequence 
logp^ coincides with the a.e. limit, i.e. with logp. By virtue of (|5.28p . we can invoke the Severini-Egorov 
theorem and deduce that for all e > there exists a set C Q such that niea.s{Q^) < e and 

logpT- logp uniformly in Q\Qe. (5.29) 



< C, 
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In addition, since (|5.6p holds, a standard interpolation calculus and the continuous embedding V C L^{ft) 
yield 



|logPr|i^{/3(Q) < 



L3(a) 



< l|logPr|li''i(o,T,Li(0)) l|logP^llL2(0,T,i'=(n)) 

< ||logPr|li''i(o,T,Li(n)) l|logPr|li2(o,T,V) ^ C, 



whence 



lll0gPr|lL«/3(Q) < C. (5.30) 

By virtue of (|5.29p - (|5.30p and of the Holder inequality, ifl<g<8/3we can write 

/ llogPr - logpl"^ < ( / |l0gp^ - logpl*^^^ ) £^~^ + / |logp^ - logpl' , 

JQ \JQe J JQ\Qe 

and then 

3£ 

hm / llogpT — logpl' < sup ( / |logpT- — logpl^^"^ ) e^~"^ < Ce^^"^. 
Taking the limit as e \ we obtain logpr — > logp in L'^{Q) VI < 9 < I, and in particular 



3' 

'2/ 



logp.^logp inL^(Q), (5.31) 

as r \ 0. Then, owing to Proposition 17.11 the same convergence easily holds (see p.2p and p.4p l for 
7^1og(p^). Besides, we can identify the weak limit y (see (15.221 1: the convergence in (|5.3ip definitely 
yields y = log p. 

Now, by the above convergences, we are in the position of taking the limit as r \ in (|4.1l) - (l4.6p and 
(1131), ((OII)) - ((iTT|l . At first, we observe that (see ([05)) and ((X5|) ) 

Pr{0) = Pqt Po in H. (5.32) 

Thanks to (g^l), (|iTn)) - (|iTT|) . ([O^ . (j^TS)) and (I^TU)) . the limit functions x, e, u satisfy the initial 
conditions (|1.4ip . (|2.3ip - (|2.32p . The rest of the proof will proceed in five steps. As first, we will pass to 
the limit in (j4.3p ; secondly we will prove that the following strong convergence holds 

dtXr-^Xt inL2(o,T,y). (5.33) 

Then, by virtue of (|5.33p we will discuss the passage to the limit in (|4.ip . Finally we will consider (|4.5p 
and we will easily recover (|2.28p by virtue of the above convergences. 

In order to pass to the limit in (|4.3p . we first observe that (j5.26p combined with (|2.39p yields (see 
Proposition 17. ip 

h{rr9r) h{9) in L'^iQ). (5.34) 
Thus, by virtue of (|5.8p . (|5.17p . (|5.3ip and (|5.34p . in order to conclude the passage to the limit as r \ 
in (|4.3p . it suffices to control the sequence ^r- In particular we need to verify that converges, in a 
suitable sense, to some selection ^ g dv,v'J{x)- To this aim, we observe that, by a comparison in ()4.3p . 
([2:391) . (EH) and dSH) imply 

lk'^llL2(o,T,y') - ^' (5.35) 

and consequently we have 

^r^^ mL\0,T,V'). (5.36) 
Hence, by (|5.17p and (|5.36p we can deduce 

(?"r,Xr>-> / (^,X>, (5.37) 



as T \ 0, which enables us to apply the result presented in [5} Prop. 2.5, p. 27] for X = L^{Q, T, V) and 
deduce 

Cedv.v'Jix) a.e. in(0,r). (5.38) 
Indeed, it is known that a maximal monotone operator from V to V' induces an analogous operator from 
L^(0, T, V) to L^(0, T, V) which is defined by the a.e. relation in (0, T). Moreover, in our framework we 
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can refer to [51 Ex. 2.3.3, p. 25] and deduce that the mduced operator from L'^{0,T,V) to L'^{0,T,V') 
is maximal monotone as weU. Finally, by the above arguments we can pass to the limit in (|4.3|) and get 
(|2.26p solved by x, p in V and a.e. in (0,T). Now. we aim to prove (|5.33p . We first note that, since 
(|5.8p holds we can infer that 

dtXr^Xt mL\0,T,V). (5.39) 
As a consequence, (j5.33p can be obtained just by verifying that 

limsup ||aar|li2(o,T,y) < IIXtlli^jo^T.v) ' (5-40) 

r\,0 

since the strong convergence of the norms combined with the weak convergence imply the required strong 
convergence (|5.33p . To obtain (|5.40p we test (|4.3p by dtXr, integrate over (0,T), and take the limsup as 
T \ 0. We have 



limsup / ||9tXr||H + ll^^tXrlln = limsup - / VXrV(9tXr - / {^.r^dtXr) 
r\,0 Jo t\0 y JQ JQ 

+ f h{Trer)dtXr~ f (logp^^tXr V (5.41) 
JQ JQ / 

We first note that, since ((5^ . (|5.17|) . (|5.34p and (|5.3ip hold, we can infer that 

lim (- / VxrVatXr + / h{TrOr)dtXT - / 7;(logp^)9tXr ) 
^N" V Jq Jq Jq J 

= - I VxVxt + / Ke)xt - I (logp)xt. (5.42) 
jq jq jq 

Next, we have to treat the term — /o^(Ct, ^tXr)- To this aim, we observe that by definition of subdiffer- 
ential and owing to (|3.ip - (|3.2p . we can write 



N N 



I {ir,dtXr) = J2i^\x'-X'-') > E./(X^) - Jix'-') = JixAT)) J(xo) (5.43) 

for any r > 0. Hence, by the lower semicontinuity in V of the function J and due to (j5.12p . we claim 
that 

limsup-/ i^r.dtXr) <~JixiT)) + J ixo). (5.44) 



r\0 

Hence, by combining (|5.42p and (|5.44p we get 



limsup - / VXrV^tXr - / i^T^dtXr) + / HTrOr)dtXT - / T; (logPr )9tXT 

r\o \ Jq Jo Jq Jq 

< - I VxVxt + / h{e)xt - I (logp)xt - J{x{T)) + J(xo), (5.45) 
Jq Jq Jq 

and the right hand side of (|5.45p is equal to 

WxtWl + vW'^XtWl, (5.46) 

as one can easily verify by testing (|2.26p by xt and then integrating in time. Thus (|5.33p is proved. 
Remark 5.1. The last result follows once proved that 

- J(X(T)) + J(xo) = - / {tXt). (5.47) 



This can be obtained by extending the statement in [51 Lemma 3.3, p. 73] to the case of abstract subdif- 
ferential operators defined in the duality pairing between V' and V . 

Now, we discuss how to perform the passage to the limit in (|4.ip . By virtue of (|5.19p we have 

Adr^Ae in L2(0,T,F')- (5-48) 
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In addition, owing to (|5.34p and (|5.33p we may infer that 

h{Trer)dtXr ^ Hojxt m l\o,t, L^^\n)) , 

dtxl~^x\ mLHO,T,H). (5.49) 
Thus, by a comparison in (|4.ip we can deduce that 

dtcr^f] in i^(0,T,F'), (5-50) 

where 77 = —AO — h{0)xt + ■ In order to show that rj = et we argue as follows. We consider the time 
convolution product 1 * (dter) which satisfies 

1 * (dtCr) = Br - eo. (5.51) 

Then, we observe that the following convergences hold thanks to (|5.50p 

1 * {dter) ^ I * V in W^i'i(0,r, V), 

while 

er-ea^e-eo in L°°(0, T, iJ). 
The uniqueness of the limit of ()5.5ip entails 

6 = 60 + 1*77. (5.52) 

Hence, e must be derivable and, by deriving (|5.52p . we get et = rj. Thus, we can pass to the limit in (|4.ip 
getting ([224]). 

Finally, we aim to take the limit in (|4.5p - (|4.6p . On account of (|5.2ip it is easily verified that 

Bpr^Bp in L°°{0,T,V')nL^{0,T,H). (5.53) 

On the other hand, thanks to (|5.10p . we have 

dtur^ut in L'^{0,T,H). (5.54) 

As a consequence of (|5.53p and (|5.54p we can pass to the limit in (|4.5p and get (|2.28p . Finally, we 
equivalently rewrite (|4.6p as follows 



1+Xr " \l+Xr 

Let us note that (|5.17p . (|5.2ip and the fact that Xr G L°°{Q) allow us to take the limit as r \ in (|5.55p 
and eventually get (|2.29p . Then Theorem 12.31 is completely proved. 

6. Regularity results 



This section is devoted to the proof of Theorems l2.4[ I2.5l and [2y6l The proof of the improved regularities 
is based on formal estimates performed on the solutions of the complete dissipative model (|2.24p - (|2.29p . 
Actually, the following estimates can be made rigorous and, in particular, the proof of Theorem 12.51 can 
be directly reproduced on the discrete scheme. 

As first, we prove Theorem 12.41 that states a further regularity for the inverse of the pressure p~^. We 
first achieve an improved regularity for the variable by means of a formal estimate, then we will 
deduce the required result for p~^. 

Before proceeding, for n < 3, we recall the Sobolev embedding H^{il) ^ L^(fi) and the Gagliardo- 
Niremberg inequality (cf. |21|). yielding, for n = 3, 

1^11^3(0) < Cgn \\v\\h llVt'll^ + Ml ■ (6-1) 

Now, we consider (|2.28p rewritten in the terms of the variables u and x as 

ut + i3(u(l + x)) -0. (6.2) 

We formally test (|6.2p by —u~^ and integrate over (0,i). For a rigorous estimate, we could proceed 
as in the First a priori estimate of Chapter 4. In particular, we could truncate the function —u~^ at 
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level e and then let e tend to zero. Let us note that, by virtue of (|2.32p . (|2.34p and (|2.52p we get 
^0 ^ = (1 + Xo)Po ^ e H. We first have 



UtU ^ 

lo Jn 

Then, by definition of B, we write 



1 



d , _2\ 1 II _i/,s||2 1 II _i||2 

^(^ )^7^h mH^olro \\h- 



(6.3) 



{B{u{l + x)),u- 







{i + x)\yu\^u 



2 „,-4 



Jn 



Jr 



(6.4) 



Jil. 



The third integral on the right hand side of (|6.4p is estimated as follows (see (|6.ip ) 



3/ f \u-'VuVx\<i f l|Vx|li6(f2)||V(u-i)||^ 
Jo Jn Jo 



lL3(a) 



< 7 Pi^~')\\Uo.t,H) + 3C^GN / llVxllie(o) \\u-'\\h ||V(?/-^)|U 





+ 3C^N / llVxIl 



(6.5) 



Now, since Vx G L°°(0, T, F) (of. p.47p ). with the help of Young's inequality, we can eventually write 



|u"^VwVx| < Cii + - ||V(w~^)| 



(6.6) 



As the trace operator 7 : — > L^(r) is compact, we may deduce that, for any > 0, there exists Cg- > 
such that 

\h{v)\\l,^^)<a\\v\\l + C,\\v\\l yveV. (6.7) 
Thus, we can control the boundary integral in (|6.4p as follows 



u-^{l + x) 



Jr 



1 1 ^ 111 II — 1 1 1 



(6.8) 



Jn 



(l + x)|Vu|'«-4 = 3 / / (l + x)|V(w-i)| 
Jo Jn 



Combining (|6.2p - (|6.8p for a sufficiently small cr, we have 



Finally we point out that for the first integral on the right hand side of ()6.4p the following equality holds 

(6.9) 

(6.10) 
(6.11) 



\u-'it)\\l + \\\7iu~')\ 



L^{0,tM} 



< Ci3 ( 1 



We can apply the Gronwall lemma [21 Theorem 2.1] to (|6.10p and finally obtain 



Ir llL°°(0,T,H)nL2(0,T,y) - ^■ 

The same estimate can be derived for the inverse of the pressure p~^. Indeed, thanks to (|2.47p and to 
the relation p^^ ~ u~^{l + x)~^, we can easily deduce that (|6.1ip implies 

11^* \\L°^{0,T,H)nL^{0,T,V) — ^- (6-12) 

Let us detail such a procedure. By virtue of the Holder inequality and of the continuous embedding 
V C L*{n), we have 



\L^^{o.T.y) 



.-1||2 

\l^{o.t,v) 



f 




2 


f 




u 


-iVx 


2 


JO 


1 + X 


H 


Jo 


1 + X 


(1 




H 



livxil 



< C 



15(11'" NL2(o_T,V) + ll^ll-L'°(o,T,iy) II" ^\\l^[0,T,V)) 



< c. 



(6.13) 
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lL°°(o T //) - II ^ IIl~(o T ff) ■ From (|6.13p we can eventually 
infer that (|6.12p holds and Theorem 12.41 is consequently proved. 



In order to prove Theorem l2.5l we perform a further a priori estimate, that can be suitably reproduced 
on the discrete scheme. Let us deal with p.24|) and rewrite it as 

(1 + xdh"{9))et + A9^ -eh\e)xt + xl 

Now, we test (|6.14p by 6t, integrate over (0,t) and get 



9l{i + xeh"m 



10 Jn 

Owing to ([2^ - ([2^ we obtain 



lrt|l2^2(o,t,_f/; 



< 



{i + x0h"ie)) 



oh'{e)0txt 



Jn 



Jn 



(6.14) 
(6.15) 



(6.16) 



Jn 



and 



Jn 
<Cie 

c, 

< - 



XtOt 



Jn 



<Ch 



\Xt\\Ot\ + 



Jn 



oh'{e)etxt 
(WxtWH + Wxtfv) IMh 

+ WxtWl + CisJ^ WxtWlWxtl 



\Xt\ 



Jn 



2 

V ' 



(6.17) 



where we have exploited the Young and Holder inequalities as well as the continuous embedding V C 
L^(il). By integrating by parts in time, thanks to (|2.33p it follows that 

(6.18) 



\Mt.^o,t.H) + l^\\'^m\\H<Cl9 + C^, / WxtWH + C^s / llXtllyllXtlly, 



where for the moment we just know that ||xt||y G L^{0,T) (see (|2.47p '). However, we are going to combine 
estimate (|6.18p with another estimate which will give us more information on ||xi||y- Indeed, we consider 

1 



(|2.26p and differentiate it with respect to time thus obtaining 

xtt + vAxtt + Axt + - h'{e)et 



-pt- 



p 



(6.19) 



Then, we test (|6.19p by xt- After some integrations by parts in time and owing to (|2.39p . we can write 



1 1 (ixtl' + HVxtl') {t) + 



|Vxt 



Jn 



<C2o||Xt(0)!|y+C„ 







1 


/ \Ot\\xt\^ 






In 


Jo Jn 


P 



(6,Xt) 

\Pt \ Ixtl 



(6.20) 



Let us note that since xo G D{dv,v' J), there exists G dv,v' J{xo)- Thus, we can introduce x'o the 
initial value of the time derivative of x by defining it as the solution of the following elliptic equation 

Xo + lyAx'o = --^Xo - & + h{9o) - logpo, (6.21) 
where the right hand side is known in V' thanks to (|2.33p - (|2.35p . (|2.39p and (|2.52p . Then, it easily follows 
that x'o € V. Next, we observe that monotonicity arguments yield 

''(6,Xt)>0. (6.22) 



By applying Holder's and Young's inequalities, we can estimate the first integral on the right hand side 
of ((5^ as follows 



Ch 



IXtllff 



(6.23) 



Analogously, the last integral on the right hand side of (|6.20p can be handled as follows 

\\Pt\\H\\Xt\\v: (6.24) 



t 

Jn 



1 




1 




\Pt\\Xt\< f 




p 


Jo 


P 



L'i{n) 
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where, owing to Theorems [13] and (||l/p!|j,4(n) ||pf !|^) G ^^0,^). Thanks to (|Og)l 
yields 



^ £(lx*l' + HVx*|') W+f /j^^'l 







1 




/' 








p 



lb*ll^llx*lly. 



(6.25) 



FinaUy, by combining (|6.18p and (|6.25p . we obtain 

\MUo,t,H) + ^ \Mml + {\xt\' + y \^xtf) (t) + II Vxt|li.(o,,,^) 



4 



<C23 + c,s \\xt\\v\\xt\\v + C24 ||x*IIh + 





1 


Jo 


p 



\\Pt\\H WXtWv 



(6.26) 



Now, we can apply the Gronwall lemma [H Theorem 2.1] and owing to Theorem 12.31 (see (|2.45p and 
dHZI)) we get 

ll^llffi(o,T,_f/)nL~(o,T,v) + \\x\\w^.=°{o,T,v) - (6.27) 

so that Theorem 12.51 is completely proved. 

Finally, we aim to prove Theorem 12.61 i.e. to establish the positivity of the temperature. To this 
purpose, we deal with (|6.14p and formally test it by —9^^. After integrating over {0,t), for the first term 
we get 



/ 0,0-'=- f f |(log0) 
Jn Jo Jn 

{logent)- f {\og9)+{t)+ I log 00, 



(6.28) 

where (log6')~ and (log 6*)+ denote the negative and positive parts of the function log 6*, respectively. 
Hence, we have 



|log0|(t) 



h"{e)0tx + 



Jn 

2 / {\og9) + {t)- I log 00 + 



Jn 



Vlog( 
h'{9)xf 



Jn 



e 



(6.29) 



Jn 



Now, we deal with the second integral on the left hand side of (|6.29p and, after integrating by parts in 
time, we can write 



h"ie)etx 



Jn 



t I- Q 

{h\e))x 



Jo 



dV 



= - I h\e{t))x{t)+ I h'{0o)xo+ / / h'{0)xt, 

Jo Jn 



(6.30) 



and then 



|log0|(t) 



|Vlog( 



Jn 



Xt 



Jn 



= / h'{e{t))x{t) + 2 / {\oge)+{t)- / log^o- / h'{e„)xo- 

Jn Jn Jn Jn 

Let us observe that well-known properties of the logarithm function and (|2.45p yield 

\^og0{t)\ 

< sup / |0|'< ||0||i=o(o,T,//) <C. 

o<t<T J nn{e>i} 



(6.31) 



(log6i)+ , 1, < sup / 

(0,1, L (u)) o<t<TJnn{9>i} 



(6.32) 
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Besides, we note that the first integral in the right hand side of (|6.3ip is bounded thanks to (|2.39p and 
((Z47| . By virtue of ((04| . ((2?57|) and ((02)l we can eventually infer that 

2 



l|logeWllLi(Jl) + l|Vlog0||i.(o^j.^^) 



Xt 



V0 



< c. 



(6.33) 



As a consequence of (|6.32p . (|6.33p and exploiting the Poincare-Wirtinger inequality, we can deduce 

2 



l|log^llL°o(0,r,Li(S^))nL2(0,T,V) + 



Xt 



< c. 



(6.34) 



L^(0,T,H) 



7. Appendix 

In this appendix, we present two auxiliary results among those exploited in our proofs. As first, 
let us detail and prove a convergence result we have applied in order to pass to the limit in the time 
discretization scheme. 



Proposition 7.1. Let Tt be the translation operator defined in 

Vr^v m L'^{-T,T,H), 
as T \ 0, then the following convergence 



If the sequence {vr} satisfies 



\\TrV,{t)~v{t)\\ijdt~~^0 



(7.1) 
(7.2) 



holds. 



Proof. Here, we give just an outline of the proof. A first step consists in proving that the following 
convergence 



\\TrW{t) - w{t)\\irdt > 

t\0 



(7.3) 



holds for any w G C°([-r, T],H^ To this aim, let us note that, if w e C°{[-T, T],H), then easily 
follows by suitably applying the Lebesgue's dominated convergence theorem. Now, since C'^{[—T,T],H) 
is dense in L^(—T,T,H), if w S L'^{—T,T,H), there exists an approximating sequence C 
C^{[-T,T],H) such that 

(7.4) 
(7.5) 

(7.6) 



— > V in L''{—T,T,H) 
as n — > +00. In particular, for any e > 0, there exists G N such that the inequality 

IK'ji — ''^\\l'^{-T,T,H) — ^ 

holds for all ji > ji^. Moreover, thanks to (|3.4p . from (|7.5p we can easily deduce that 



/ WTrVnit) - TTv{t)\\\dt < ||w„ - v\\l2^_TT H) ^ 

Jo 



for all n > and for all < r < T. Now, we fix n = and point out that since Vn G C°([— T, T], H) 
then (|7.3p holds for w = Vn- In particular, there exists such that 

WTrVnit) - Vnit)\\l dt < e (7.7) 



for all T < Tg. This allows us to split the difference (TtV — v) into three terms, all tending to zero in 
L^{0,T, H). Namely, we can write 



TtV — U = TtV — TtV„ 
TtVji Vji 
+ Vn-V. 



(7i 



^By C^{[—T,T],H) we denote the space of continuous functions from [— T,T] into H equipped with the L°° {~T,T, H) 
norm. 
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By means of the above arguments, we can infer that the foUowing inequahty 

r-T 







\Trv{t)-v{t)\\]jdt<ie (7.9) 



holds for all r < r^, whence 



\\Trv{i)-v{t)\\ldt~^Q. (7.10) 
Finally, we consider the difference between %-Vr and v. It is convenient to express it as follows 

TtVt — V ~ TtVt ~ TtV 

+ TrV~-V. (7.11) 

Thanks to (|7.ip and to the definition of for the first addend on the right hand side we can write 

WTrVrit) - TrV{t)\\% < \\Vr - ^11 i2(_T,T,ff) ^^^^ 0' (7-12) 

for all < T < T. While the second addend has already been discussed in ()7.10p . Hence, we can infer 
that 

\\TrVr{t)-v{t)\\l^dt >0 (7.13) 

and the proof is complete. □ 

Finally, we present a fairly standard result concerning the sum of maximal monotone operators. Let 
us recall some notations previously introduced. We denote by H the space L^(ri), with a bounded 
domain included in M.^, and by A the abstract operator prescribed in (|2.10p . Next, by j we denote a 
proper, lower semicontinuous and convex function on H, and by dnj the subdifferential oi j. Then the 
following result holds. 

Theorem 7.2. Let f G H. Then, there exists a unique x G H'^{Vl) n D{j) such that 

X + Ax + dH]{x)^ .!■ (7.14) 

The above result can be easily proved by approximating the subdifferential operator in (j7.14p by its 
Yosida regularization {dHj)x- We remind that {dHj)x is maximal monotone and Lipschitz continuous. 
Hence, we may invoke well known results on the sum of maximal monotone operators (cf. e.g. |8]) and 
get the existence of a solution x\ € H^{fl) solving the approximated equation. Next, we could perform 
standard a priori estimates, independent of A, on the approximated system and get the limit as A \ by 
compactness. Finally, uniqueness could be deduced via standard contradiction arguments owing to the 
monotonicity of the subdifferential operator. 
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